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Effects of geometrical frustration in low-dimensional charge ordering systems are theoreti- 
cally studied, mainly focusing on dynamical properties. We treat extended Hubbard models at 
quarter-filling, where the frustration arises from competing charge ordered patterns favored by 
different intersite Coulomb interactions, which are effective models for various charge transfer- 
type molecular conductors and transition metal oxides. Two different lattice structures are 
considered: (a) one-dimensional chain with intersite Coulomb interaction of nearest neighbor 
Vi and that of next-nearest neighbor V2, and (b) two-dimensional square lattice with Vi along 
the squares and V2 along one of the diagonals. From previous studies, charge ordered insulat- 
ing states are known to be unstable in the frustrated region, i.e., Vi ~ 2V2 for case (a) and 
Vi ~ V2 for case (b), resulting in a robust metallic phase even when the interaction strenghs are 
strong. By applying the Lanczos exact diagonalization to finite-size clusters, we have found that 
fluctuations of different charge order patterns exist in the frustration-induced metallic phase, 
showing up as characteristic low energy modes in dynamical correlation functions. Comparison 
of such features between the two models are discussed, whose difference will be ascribed to the 
dimensionality effect. We also point out incommensurate correlation in the charge sector due 
to the frustration, found in one-dimensional clusters. 

KEYWORDS: charge order, geometrical frustration, molecular conductors, transition metal oxides, exact 
diagonalization, extended Hubbard model 



1. Introduction 

Charge ordering (CO) phenomena have attracted in- 
terest from early days, motivated by experiments on vari- 
ous materials. A classic example is transition metal oxide 
Fe304, the magnetite, where its metal-insulator transi- 
tion was proposed by Verwey to be due to CO among the 
electrons of the Fe ions at the so-called i?-sites.^'^ How- 
ever, this issue has not been solved despite enormous 
amount of works, and even the existence of CO itself 
has recently been doubted. '^^^ A reason for such puzzle 
is that this CO system is under strong geometrical frus- 
tration as pointed out by Anderson.^ This is seen when 
the system is described by Ising variables cr = ±1 corre- 
sponding to occupied/unoccupied sites, on the pyroclore 
lattice formed by the B-sites. An antiferromagnetic in- 
teraction for each nearest neighbor (ij) bond arises due 
to the Coulomb repulsion favoring aiaj — —iJ This is 
indeed a typical situation of geometrical frustration, as 
has frequently been discussed in spin systems.* 

CO systems under geometrical frustration have re- 
newed its interest rather recently, triggered by differ- 
ent experimental works. The one-dimensional (ID) Cu-0 
unit in a transition metal oxide PrBa2Cu408^ has been 
pointed out to be susceptible to frustration due to its 
zigzag chain structure.^'' A class of molecular conductors, 
i.e., e-ET2X, where FT is an abbreviation for BEDT- 
TTF (bisethilenodithio-tetrathiofulvalene) molecule and 
X stands for different monovalent anions with closed 
shell, "'^^ is a clear example of systems affected by frustra- 



tion in their two-dimensional (2D) layer. Similar situ- 
ation is realized in many molecular conductors as well, 
generally described by 2D anisotropic triangular lattice 
models at quarter-filling. ^'^"^^ Such interplay between 
CO instability and geometrical frustration has also been 
discussed in triangular lattice systems such as LuFe204 



18 



and Naa;Co02, spinel compounds such as AIV2O4 
and LiV204,^^ and so on. 

In these compounds, the CO transition, when it is 
realized, takes place accompanied with a large mod- 
ulation in the lattice degree of freedom. For exam- 
ple, a first order structural phase transition is observed 
in 6'-FT2RbZn(SCN)4, where the ET molecules show 
displacements from their positions at high tempera- 
tures.^^' Ordering of Na"*" cations is suggested to be 
closely related to the CO state in Nao.5Co02.^^' In 
AIV2O4, a pecuUar type of CO with a 3:1 ratio of Y'^-^+s 
and V^ '^^'^'^ ions has been suggested from a large rhom- 
bohedral distortion.^'* These instabilities can be consid- 
ered as "efforts" of the system to relax the geometrical 
frustration and to settle down in stable ordered states. 

On the other hand, when the compounds avoid such 
structural instabilities, conducting states are realized. 
PrBa2Cu408 is metallic down to the lowest tempera- 
ture^ but large charge fluctuation is seen in the NQR 
relaxations.^'' In 6'-ET2CsZn(SCN)4 and in the rapid- 
cooled 6'-ET2RbZn(SCN)4, the electrical conductivity 
shows rather little temperature dependence down to 
about 20 K, below which a gradual increase of resitivity 
is observed.^-' There, an NMR measurement^^ suggests 
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the existence of a "charge glass" -Hke state where the 
dynamical CO fluctuation becomes frozen at low tem- 
peratures but not in a long-range way. Moreover, non- 
linear conductivity and insulator-like dielectric behavior 
are observed even when the systems are in the conductive 
state.26 

In this paper, we theoretically investigate such CO 
systems under strong geometrical frustration, by consid- 
ering models which have been discussed previously to 
describe some of the materials above, i.e., geometrically 
frustrated quarter-filled extended Hubbard models which 
include competing intersite Coulomb interaction terms 
Vij in addition to the on-site [/-term. Previous numerical 
works^°'^^'^^'^* have revealed that CO insulating states 
in the unfrustrated region are melted when the frustra- 
tion becomes strong, resulting in a metal even when the 
interaction strengths, U and Vij, are large. In this metal- 
lic state with large but frustrated Vij, qualitative differ- 
ence has hardly been found, mainly in static properties of 
the models, in comparison to a metal at "weak-coupling" 
regime with small Vij (but with large U). We will see 
in the present paper that differences exist in dynamical 
quantities. Charge fluctuations of different CO patterns 
exist at finite energy in this metallic region, which is a 
novel property of CO system under strong geometrical 
frustration. 

The organization of this paper is as follows. In § 2, 
ID and 2D frustrated extended Hubbard models which 
we consider, a summary of previous works on them, and 
the numerical Lanczos exact diagonalization method we 
use in this paper, are introduced. The results on the ID 
and the 2D cases are shown in § 3 and 4, respectively, 
in which mainly the dynamical correlation functions are 
discussed. Section 5 is devoted to discussions and § 6 to 
a summary. 

2. Model and Method 

We consider quarter-filled extended Hubbard models 
on two different lattice structures. The first is a ID chain 
with nearest neighbor and next-nearest neighbor interac- 
tions, described as, 

Hid '^{-'tl{cl,cr'^i+i,<^ + h.C.) - t2(c|,o-Ci+2,o- + h.c.)} 

+ ^ {Un^^irii^i + ViTiiTii+i + V2nini+2) , (2.1) 

i 

where cr is a spin index taking "f and |, {ci,a) and 

ni^a = cl ^Ci^a denote the creation (annihilation) opera- 
tor and the number operator for the electron of spin a 
at the i-th site, and rii ~ ni,^ + rii^i. The other is a 2D 
square lattice model with next-nearest neighbor interac- 
tion along one of the diagonals, 

H2D = ^{-ii(4^Ci+£,CT -I- cl^Ci+y^o- + h.c.) 

ia 

- t2{c\,^Ci+x+yM + h.c.)} 

+ ^ {C^".i,T"*,i + ^1 + niUi+y) + V2^^^«^+^+^l} i 

(2.2) 
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Fig. 1. (Color online) (a) One-dimensional chain with next- 
nearest neighbor interaction (up) and zigzag ladder (down), (b) 
square lattice with next-nearest neighbor interaction along one 
diagonal (left) and anisotropic triangular lattice (right). These 
models, respectively, are equivalent with each other as seen in 
the figure. 

where x and y are the unit vectors with one lattice spac- 
ing along X and y directions, respectively. As shown in 
Fig. 1, these models are equivalent to a ID double-chain 
model with zigzag couplings between the chains, i.e., the 
zigzag ladder, and a 2D model on an anisotropic trian- 
gular lattice, respectively. The lattice constants in both 
models are set to unity in the following. 

The former ID model at quarter-filling has been stud- 
ied as an effective model for two different systems. One 
is for quarter-filled molecular conductors TMTSF2X, 
where effects of long-range Coulomb interaction have 
been investigated.^^' '^'^ In these studies ti is set to while 
V2 < V\ is appropriate to the actual compounds since 
the Coulomb interactions are taken along the TMTSF 
chain. The other is for the Cu-O subunit of PrBa2Cu408, 
mentioned in § 1, where the Cu sites are formed in a 
zigzag ladder structure. In this case \t\\ ^ \t2\ due 
to the nearly 180° (ii) and 90° (^2) Cu-O-Cu angles, and 
V\ ~ \/2V2 suggested by the Cu-Cu distances. For ti = 
or t2 = 0, this ID model at [/ = 00 is mapped onto XXZ 
models^° therefore the analogy between CO systems and 
localized spin systems under frustration metioned in § 1 
is straightforward here. 

The latter 2D model, on the other hand, together 
with its variants with lower symmetry, has been stud- 
ied intensively motivated by recent experimental find- 
ings of CO in 2D molecular conductors. ^^'^'^'^^'^^'■^-^'•^^ 
Among them, compounds with the highest symmetry, Q- 
ET2Ar, can be described by the model in eq. (2.2). Here, 
X takes different monovalent anions, e.g., MM'(SCN)4 
with various metal elements in MM' such as M=Rb, Cs, 
Tl and M'=Zn, Co." In these 6'-ET2MM'(SCN)4 salts, 
the anisotropic shape of the HOMO forming the valence 
band at the Fermi level results in ti ^ t2 while relative 
distances between molecules provide Vi ~ V2. 

Numerical studies^°' ^^'^^'^^ on these ID and 2D mod- 
els have found a wide region of metallic phase between 
two CO insulating phases in their ground state phase 
diagrams, as shown in Fig. 2. This is drawn by varying 
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Fig. 2. (Color online) Schematic ground state phase diagram of 
frustrated extended Hubbard models, for fixed on-site Coulomb 
energy [/ at a large value compared to the kinetic energy tij. 
COl and C02 denote the charge ordered insulating phases fa- 
vored by the competing intersite Coulomb repulsions V\ and V2 1 
respectively, and M is the metallic phase. (yi)c and (V2)c are 
their critical values for the non-frustrated cases. Note that the 
detailed nature of the phases and their boundaries can be differ- 
ent depending on the lattice structures; for example the charge 
ordering and the metal-insulator transitions take place at differ- 
ent parameters in the 2D models.^*' '^'^ 



Vi and ( < U /2) while U is fixed at a large value 
compared to the kinetic term tij. In the two CO phases, 
Wigner crystal-type CO states are stabilized due to ei- 
ther Vi or V2 : charge localization occurs in order to avoid 
the Coulomb repulsion along the nearest neighbor or the 
next-nearest neighbor bonds. In the ID case, they are 
the "[1010]" and "[1100]" (or 4fcF and 2A:f CO states in 
terms of fcp — 7r/4 for <2 = 0) patterns, while in the 
2D case, the checkerboard-type and the stripe-type pat- 
terns, respectively. The metallic phase is realized in the 
Vi ~ aV2 region, where a = 1 (a = 2) for the ID (2D) 
model, due to strong geometrical frustration in Vij . The 
line along Vi = aV2 is the "fully" frustrated line where 
the corresponding Ising models Vijataj have a spin 

disordered ground state due to macroscopic degeneracy. 

In this work, in order to incorporate quantum fluc- 
tuation which is crucial in describing this frustration- 
induced metallic state, we use standard numerical Lanc- 
zos exact diagonalization method for finite size clus- 
ters.'^'* The electron number n is fixed at quarter-filling, 
namely, n — L/2 where L is the total number of 
sites in the cluster. We mainly consider the dynamical 
charge/spin correlation functions, 

N{q,Lo) = Y,\{^\Ng\0)\^S{LU -E, + Eo) (2.3) 

S{q, ('^I^Q |0) -^- + ^0) (2.4) 

where i?o and E^j are the ground state (|0)) and excited 
state {\v)') energies of the system, respectively. TVg and 
are Fourier transformations of the electron density 
with respect to the mean value h — n/L = 1/2, Ni — 
Hi — n, and the z-compenent of the spin operator, — 
Nq = and S*^ = 



{c\^Ci^-c\^Cii)/2, I.e. 
In order to consider the CO patterns expected from 



the previous works, including the possible spin patterns, 
L = 4m (to: integer) site chains are required for the 
ID case, while in the 2D case the L = 4 x 4 = 16 
site square-shaped cluster is the only available size when 
considering the matching of the cluster shape with the 
periodicities of such CO/spin order, and the computa- 
tional limitation; for the ID case the L — 16 site chain is 
used. The antiperiodic boundary condition is choosen for 
both cases to give the closed shell in their non-interacting 
bands. '^^ Note that the latter 2D cluster rules out long- 
period charge modulations recently found in the frus- 
trated metalhc state in this 2D model.^^'^^'^^ 

In the following calculations we set ti as unity, and fix 
t2 = and U = 10 while varying the intersite Coulomb 
energies Vi and ¥2- Although some of the properties of 
the models might be modified by different values of t2, 
as far as the properties discussed in this paper we expect 
that they will not distract our main conclusions. 

3. One-dimensional case 

Before discussing the dynamical properties of the ID 
model, let us remark a result on static properties that was 
not noticed in the previous studies: an incommensurate 
(IC) peak in the equal-time charge correlation function, 
N{q) = i-i(0|7V_q7Vq|0) = L-^J:^Je"'■''Hm^NJ\0), 
where Rij denotes the vector connecting Rt and Rj, 
when the system is in the frustration-induced metallic 
phase. In Fig. 3, profiles of N{q) for fixed Vi = 4 are 
plotted, for different values of V2. For this value of Vi, 
by increasing the value of V2 from 0, the ground state 
changes from the [1010]-type Akp CO insulator to the 
metal in the intermediate region, and then to the [1100]- 
type 2/eF CO insulator. These two CO states are char- 
acterized by commensurate N{q) peaks at q = tt and 
7r/2, typically seen in Fig. 3 for V2 = 1 and V2 = 4, re- 
spectively. However, for V2 = 2 one can see a double peak 
structure at q = tt and at an IC value q ~ 37r/4 (which 
may be a convolution of q = tt and q ~ 57r/8 peaks). For 
V2 = i the IC peak shifts to q ~ Stt/S and dominates 
over the q = n peak. 

We note that this result does not immediately point 
to the existence of an IC charge density modulation in 
the ground state of this model. Although we find such 
IC peaks also in different cluster sizes, our study is un- 
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Fig. 3. (Color online) Profile of N{q) in the one-dimensional frus- 
trated extended Ifubbard model (L = 16) at quarter-filling, for 
fixed U = 10, Vi = 4 and for different values of V2 . 
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ID model (L=16), U=10, V^=0 



ID model (L=16), U=10, ¥2=0 
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Fig. 4. (Color online) N{q,u!) for the quarter-filled one- 
dimensional extended Hubbard model (L = 16) without frus- 
tration (V2 = 0) at fixed U = 10, for different values of Vi. 
The system size is L = 16 and the antiperiodic boundary condi- 
tion is used. The Lorentzian broadening of the delta functions is 
7? = 0.1. 



sufficient for investigating the long distance behavior 
of this IC correlation, since considerably large system 
sizes are needed to make many q values available and 
to proceed finite size scaling of them. Considering the 
density-matrix-renormalization-group results by Nishi- 
moto et al?''^ supporting the Tomonaga-Luttinger liquid 
(TLL) properties in this metallic regime, we expect that 
this IC correlation would decrease exponentially and van- 
ish at long distances since the TLL implies that low-lying 
excitations are described by modes at q = mk-p (jn: inte- 
ger). Still, this result is notable since similar long period 
modulations are found by different authors in the 2D 
models in the frustrated region,^^"'^^''^^ although we can- 
not find them in the 2D cluster that we study as pointed 
out in § 2. We will mention about it in § 5. 

Next, we proceed to the dynamical properties of this 
model. For clarity, we first show how the quantities be- 
have in the case without frustration, and then switch on 
the frustration. In Figs. 4 and 5, N{q,uj) and S{q,uj) at 
V2 = 0, i.e., the usual ID extended Hubbard model, for 
several values of Vi are shown. For Vi = V2 = 0, the 
model is just the ID Hubbard model and the obtained 
spectra are consistent with the large U study by Ogata 
and Shiba^^ based on the Bethe ansatz, as well as those 
of the ID t-J model, -^^'^^ which can be considered as an 
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Fig. 5. (Color online) S{q,w) for the quarter-filled one- 
dimensional extended Hubbard model (L = 16) without frus- 
tration (V2 = 0) at fixed U = 10, for different values of Vi. The 
Lorentzian broadening of the delta functions is = 0.1. 



effective model for the Hubbard model at large U. The 
spectra are composed of low-lying modes in the charge 
sector at Akp and in the spin sectors at 2/cf,'^^ i-e., at 
q = n and at q = 7r/2 for quarter-filling, respectively. 
When Vi is increased, in N{q,uj) a strong softening of 
the q = TT mode is seen, signaling the instability toward 
the 4/cf CO. At the same time, in S{q,uj) the low-lying 
q — tt/2 mode becomes sharper and around q — nil 
the dispersion is more symmetric with respect to this q 
value. One can interpret this behavior as a crossover to a 
more localized spin-like state due to the charge localiza- 
tion, with a well-defined spin-wave dispersion centered 
at q = 7r/2. The difference in the energy scale for the 
charge and spin degrees of freedoms is noticeably seen in 
the dispersions of these modes as distinct total widths of 
them, which is due to the fact that the CO is determined 
by the Coulomb energies V and Y while the spins are in- 
teracting with an effective Heisenberg coupling, J.'^^ 

The variation of these spectra by the inclusion of frus- 
tration controlled by the value of V2 is shown in Figs. 6 
and 7. While fixing at = 4, by increasing V2 the charge 
fiuctuation seen in A^(q, ui) shows a drastic modification. 
As seen in Fig. 6, the spectrum for V2 = 4, at which 
the system is expected to be in the region of the 2fcF 
CO state, a low-lying mode is mostly concentrated at 
q — -K 12. This is similar to the case of the 4A:f CO state 
that we have observed in Fig. 4 with the q = tt mode. In 
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ID model (L=16), U=10, Vi=4 
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Fig. 6. (Color online) N{q,u)) for the quarter-filled one- 
dimensional extended Hubbard model (L = 16) with frustration 
at fixed U = 10, Vi = 4, and for different values of V2- The 
Lorentzian broadening of the delta functions is Jj = 0.1. 
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between, in the frustration-induced metallic region for 
V2 = 2 and for V2 = 3, we can see that the spectrum 
smoothly crosses over between the dispersions charac- 
terstic of the two different CO patterns. Interestingly 
enough, prominent modes at different q values are simul- 
taneously seen in this region. For example, for V2 = 2, 
peaks at q = TT and q = Stt/S are larger compared to 
the other neighboring q values. As for the spectrum for 
V2 — 3, the q — it/2 peak becomes comparable to the 
latter peak. These match with the IC peak in N{q) men- 
tioned above, i.e., the integration of N{q,uj) along the 
w-axis. These results suggest that fluctuations of differ- 
ent CO states co-exist as dynamical charge correlation 
in the frustration-induced metallic region. We note that 
these are consistent with the TLL state at low energy 
since the energy positions of the dispersion always show 
minima at either q = 7r/2 or tt, matching the above con- 
dition q = mkp (m: integer), even when the peak heights 
are highest at a different q. 

The spin degree of freedom is consistent with such fea- 
tures. As seen in Fig. 7, the spectra with the spin-wave 
mode centered at q = 7r/2 seen in the V2 — case for the 
Akp CO state show a smooth crossover to the spin exita- 
tion at q = TT characteristic of the 2kp CO state at large 

29,30 ^]-^g latter state the pronounced peak is seen 
at q = TT although it stays at finite energy even when V2 
is largely increased, which is due to the expected spin 



gap, as the [1100] configuration of charge localization 
produces pairs of singlet formation. In the frustration- 
induced TLL region, the S{q, lu) spectrum resembles to 
that in the unfrustrated case at = V2 = in Fig. 5, 
in contrast with their clear difference in N{q,LL!). 

4. Two-dimensional case 

In the 2D model, the crossover behavior between two 
CO insulating states is also observed in the dynamical 
properties of the system. In Figs. 8 and 9, N{q,uj) and 
S{q, uj) spectra^" are shown for different (Vi, V2) parame- 
ter sets across the phase diagram in Fig. 2 (see also Fig. 1 
of Ref. 12). When one starts from the checkerboard-type 
CO insulating state in the (large Vi, small V2) region, by 
decreasing Vi and increasing V2, the ground state trans- 
forms to the intermediate metallic and then to the stripe- 
type CO insulating phase. The two CO states are char- 
acterized by wave vectors q — (tt, it) for the checkerboard 
pattern'^'^^ and q = (tt, 0) for the stripe pattern. 
These are seen in the large peaks at ~ in the N{q, uj) 
spectra, as typically seen in Fig. 8 for {Vi,V2) = (3.5, 0.5) 
and (0.5,4.5), respectively. At other wave vectors, peaks 
are relatively small and located at higher energies. 

On the other hand, in the frustration-induced metal- 
hc phase at (Vi,V2) = (4,4), both of the two peaks at 
q = (tt, tt) and (tt, 0) are prominent at the same time, but 
located at higher energies than those in the CO states. It 
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Fig. 9. (Color online) S{q, u>) for the quarter-filled two- 
dimensional extended Hubbard model (L = 16) with frustra- 
tion at U = 10, and for different sets of {Vi,V2)- The Lorentzian 
broadening of the delta functions is ?j = 0.1. 



shows that these modes are still present, although, due 
to the hmited number of data for different wave vectors 
in contrast with the ID case, how the curve of the dis- 
persion would be in the thermodynamic limit is difficult 
to judge; actually the positions of the peaks among the 
available q values show little variation in energy. This 
is more clearly seen when we compare with the N{q, lu) 
data for = (0,0), shown in Fig. 8 for reference, 

which is more dispersive and continuum-like peaks show 
up. The clear distinction between these two spectrum in 
the metallic region leads us to conclude that charge fluc- 
tuations for different CO states co-exist in the strongly 
frustrated region. This is similar to the ID case where 
different CO modes were observed in the intermediate 
metallic phase, but with more clear dispersions (Fig. 6). 

The S{q,uj) spectra, again, are consistent with such 
crossover behavior, although its behavior is apparently 
rather complicated, fn the checkerboard CO state the 
spins are located on the charge rich sites and antiferro- 
magnetically coupled. Due to the cluster shape we take, 
the spin pattern can be described by different q values 
such as <7 = (7r/2,7r/2) or (7r,0). The latter is devel- 
opped as a pronounced low-energy peak in Fig. 9 for 
(1^1, V2) = (3.5,0.5). In our calculation it is found that 
spin pattern with either of the above q is developped 
(or both) can depend on the parameters of the model 



and boundary conditions. This is probably a finite size 
effect, but we could not specify the actual origin of it. 
On the other hand, in the stripe CO state the ID chains 
are coupled antiferromagnetically along the chains. The 
interchain configuration, this time, is in fact expected 
to depend on the parameters of the model in the ther- 
modynamic limit, suggested by previous mean- field cal- 
culations.-^^' "^^ The peaks seen in both q = (tt, 0) and 
(7r,7r) in the (Vi,V2) = (0.5,4.5) indicate that the an- 
tiferromagnetic stripes have ferromagnetic correlation in 
their transverse directions. However, again, we find that 
not only the parameters of the model but also the bound- 
ary condition can change this transverse phase, therefore 
the results about the actual q value in this small cluster 
size should not be taken conclusively. Note the differ- 
ence in the energy scale between N{q,uj) and S{q,LLj) 
spectra in these CO states, where small but noticeable 
peaks exist at higher energy in the former while almost 
all spectral weight is concentrated at w < 2 in the lat- 
ter; this is analogous to the ID case. In between, in 
the intermediate metallic region the spectrum is com- 
posed of modes surviving at the above q values, as seen 
in the {Vi,V2) = (4,4) data. The S{q,uj) spectrum at 
(Vi, V2) = (0, 0) is slightly more dispersive and extended 
to higher energy region than the rather flat dispersion 
in this (Vi,V2) = (4,4) data for the strongly frustrated 
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metal. 

5. Discussions 

The IC correlation we found in the ID model in § 3 is 
a typical consequence of geometrical frustration in gen- 
eral. We suspect that it would also appear in 2D model as 
well if one could have larger system sizes available. Such 
a state supposedly correspond to the so-called "three- 
fold" state found in mean-field studies^^'^^ as well as in 
a variational Montc-Carlo study^* of this 2D model. A 
recent RPA study^^ found that the position of q where 
the charge susceptibility shows a maximum is shifted 
by changing the range of intersite Coulomb interaction 
terms and their values, and can be adjusted to that found 
in X-ray experiments on ^-ET2X as diffusive rods, i.e., 
as a short range order. Such an IC correlation is found in 
spin systems as well. Exact diagonalization calculations 
for finite size ID clusters of the so-called Majumdar- 
Ghosh model, i.e., the 5 = 1/2 Heisenberg model with 
next-nearest neighbor coupling, show such IC peak in 
its spin structure factor S{q),'^^ and an IC spin ordered 
ground state is actually found experimentally in a quasi- 
ID S* = 1/2 spin system Cs2CuCl4.''^ It is interesting to 
search for such an IC state attaining long range order 
of charge modulation in the CO materials as well. This 
will bring about different aspects from the spin systems, 
since for the CO systems the interplay with the metal- 
insulator transition would be an interesting issue. 

We consider that the co-existence of different CO 
modes in the dynamical spectra for the strongly frus- 
trated metallic region that we have shown in both ID 
and 2D models, is a general characteristic of geometrical 
frustrated CO systems. In the 2D model, a variational 
Monte Carlo simulation by Watanabe and Ogata^^ sug- 
gested that variational states with different CO patterns, 
i.e., the checkerboard, the stripe, and the threefold- type, 
have very close ground state energies in the frustrated 
metallic region. Although their data is for static ground 
states while ours is for dynamic quantities, we consider 
that we are detecting different aspects of the same phe- 
nomena. 

The flatness of the dispersions in this region for the 
2D model shown in § 4, both in N{q,uj) and S(q,uj), is 
a peculiar point. In the ID case, in contrast, the width 
of the dispersion in N{q,uj) is not much reduced, and 
that in S{q, lu) hardly differs, compared with the "weak- 
coupling" metal at 14 = V2 = 0. This indicates that 
the charge and spin excitations are very localized in real 
space for the 2D system, while it can propagate to some 
extent in the ID case. It is due to the dimensionality 
effect, since, e.g., an addition of charge, or a flipping of 
spin in the 2D model would distract the environment 
due to the strong coupling strenghs, and therefore can- 
not propagate into the "bulk". The flatness results in a 
"downward shift of spectral weight" compared with the 
Vi = V2 = data, which is pointed out by Ramirez* to 
be a tendance of geometrically frustrated spin systems. 

The relation between our results with the experimental 
data showing the existence of charge fluctuation at very 
low energy scale detected by NQR/NMR is still obscure. 
By NQR measurements on PrBa2Cu408, Fujiyama et 
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al?^ have evaluated the fraction of such fluctuating (or 
frozen) CO to be very small compared to the bulk, and 
they have deduced this to be related with small amount 
of impurities. Kanoda et al?^ have proposed the exis- 
tence of a charge glass state in 0-ET2X, an analogy of the 
spin glass state in spin systems, therefore again points to 
the role of impurities under the presence of geometrical 
frustration. In addition, CO is in general coupled with 
the lattice, not included in our model, which may also 
be important in understanding such phenomena. How- 
ever, the behavior of the ID and 2D frustrated extended 
Hubbard models we discussed in this paper would be 
helpful in analyzing such experiments to start with. 

6. Summary 

In summary, we have theoretically investigated CO 
systems under strong geometrical frustration, by treating 
frustrated extended Hubbard models on ID and 2D lat- 
tice structures. Charge fluctuations of different types of 
CO states are found to co-exist in the frustration-induced 
metallic region of these models, showing up as dynam- 
ical modes in the correlation function spectra. We have 
also found incommensurate charge correlation in the ID 
model, and deduced that it is a general tendence in other 
frustrated CO systems as well. 
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